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We introduce a new model for dark energy in the universe in which a small cosmological constant 
is generated by ordinary electromagnetic vacuum energy. The corresponding virtual photons exist 
at all frequencies but switch from a gravitationally active phase at low frequencies to a gravita- 
tionally inactive phase at higher frequencies via a Ginzburg-Landau type of phase transition. Only 
virtual photons in the gravitationally active state contribute to the cosmological constant. A small 
vacuum energy density, consistent with astronomical observations, is naturally generated in this 
model. We propose possible laboratory tests for such a scenario based on phase synchronisation in 
superconductors. 

PACS numbers: 95.36.+x; 74.20.De; 85.25 Cp 



Current astronomical observations [E 0, i, || provide 
compelling evidence that the universe is presently in a 
phase of accelerated expansion. This accelerated expan- 
sion can be formally associated with a small positive 
cosmological constant in the Einstein field equations, or 
more generally with the existence of dark energy. The 
dark energy density consistent with the astronomical ob- 
servations is at variance with typical values predicted by 
quantum field theories. The discrepancy is of the order 
10 122 , which is the famous cosmological constant problem 
Q . A large number of theoretical models exist for dark 
energy in the universe (see e.g. [y, [7[ for reviews). It is 
fair to say that none of these models can be regarded as 
being entirely convincing, and that further observations 
and experimental tests [7|, l8|, l9( are necessary to decide 
on the nature of dark energy. 

The most recent astronomical observations Q seem 
to favor constant dark energy with an equation of state 
w = — 1 as compared to dynamically evolving models. In 
this paper we introduce a new model for constant dark 
energy in the universe which has several advantages rel- 
ative to previous models. First, the model is conceptu- 
ally simple, since it associates dark energy with ordinary 
electromagnetic vacuum energy. In that sense the new 
physics underlying this model does not require the pos- 
tulate of new exotic scalar fields such as the quintessence 
field. Rather one just deals with particles (ordinary 
virtual photons) whose existence is experimentally con- 



* Electronic address: |c.beckBqmul.ac.uk| 

URL: http: //www. maths .qmul . ac .uk/~ beck 

^Electronic address: michacl.mackcy@mcgill.ca 

URL: http://www.cnd.mcgill.ca/people_mackey.html Also: 
Institut fiir theoretische Physik, Universitat Bremen, Germany 



firmed. Secondly, the model is based on a Ginzburg- 
Landau type of phase transition for the gravitational ac- 
tivity of virtual photons which for natural choices of the 
parameters generates the correct value of the vacuum en- 
ergy density in the universe. In fact, the parameters in 
our dark energy model have a similar order of magnitude 
as those that successfully describe the physics of super- 
conductors. Finally, since the phase of the macroscopic 
wave function that describes the gravitational activity of 
the virtual photons in our model may synchronize with 
that of Cooper pairs in superconductors, there is a possi- 
bility to test this electromagnetic dark energy model by 
simple laboratory experiments. 

Recall that quantum field theory formally predicts an 
infinite vacuum energy density associated with vacuum 
fluctuations. This is in marked contrast to the observed 
small positive finite value of dark energy density pdark 
consistent with the astronomical observations. The rela- 
tion between a given vacuum energy density p vac and the 
cosmological constant A in Einstein's field equations is 

A 8nG 

A = -^Pvac, (1 

c 

where G is the gravitational constant. The small value 
of A consistent with the experimental observations is the 
well-known cosmological constant problem. Suppress- 
ing the cosmological constant using techniques from su- 
perconductivity was recently sug gest ed in a paper by 
Alexander, Mbonye, and Moffat [lj|. To construct a 
simple physically realistic model of dark energy based 
on electromagnetic vacuum fluctuations creating a small 
amount of vacuum energy density p vac — Pdark, we as- 
sume that virtual photons (or any other bosons) can ex- 
ist in two different phases: A gravitationally active phase 
where they contribute to the cosmological constant A, 
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and a gravitationally inactive phase where they do not 
contribute to A. 

Let | 'J | 2 be the number density of gravitationally active 
photons in the frequency interval [v, v + dv\ . If the dark 
energy density pdark of the universe is produced by elec- 
tromagnetic vacuum fluctuations, i.e. by the zero-point 
energy term ^hv of virtual photons (or other suitable 
bosons), then the total dark energy density is obtained by 
integrating over all frequencies weighted with the number 
density of gravitationally active photons: 



Pdark 



The standard choice of 
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|*| 2 = - . 47TI/ 2 , 
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in which the factor 2 arises from the two polarization 
states of photons, makes sense in the low- frequency re- 
gion but leads to a divergent vacuum energy density for 
v — > oo. Hence we conclude that |W| 2 must exhibit a 
different type of behavior in the high frequency region. 

In the following we construct a Ginzburg-Landau type 
theory for |\&| 2 . Our model describes a possible phase 
transition behavior for the gravitational activity of vir- 
tual photons in vacuum, which has certain analogies with 
the Ginzburg-Landau theory of superconductors (where 
l^j describes the number density of superconducting 
electrons). It is a model describing dark energy in a 
fixed reference frame (the laboratory) and is thus ideally 
suited for experiments that test for possible interactions 
between dark energy fields and Cooper pairs [§| . 

We start from a Ginzburg-Landau free energy density 
given by 



F = a|*| 
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where a and b are temperature dependent coefficients. In 
the following we use the same temperature dependence 
of the parameters a and b as in the Ginzburg-Landau 
theory of superconductivity [Til . [l2j : 



a(T) 
b(T) 
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Here t is defined as t := T/T c , T c denotes a critical 
temperature, and ciq < 0, bo > are temperature- 
independent parameters. Clearly a > 0, b > for T > T c 
and a < 0, b > for T < T c , . The case T > T c describes 
a single-well potential, and the case T < T c a double-well 
potential. 

The equilibrium state ^ eq is described by a minimum 
of the free energy density. Evaluating the conditions 
F'(^ eq ) = and F"(^ eq ) > 0, for T > T c we obtain 



0, F e , 



0. 
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whereas for T < T c 
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In the following, we suppress the index eq . 
With eqs. {5J and © we may write eqs. 

I*| 2 = -?(1-* 4 ) 



as 
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For very small temperatures (T << T c ) one has l^l 2 = 
— ao/bo, which we identify with the low- frequency behav- 
ior of photons as given by eq. ([3]) . Thus 
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which leads to 
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We also need to formally attribute a temperature T to 
the virtual photons underlying dark energy. This can be 
done as follows: Virtual photons have the same energy 
as ordinary photons in a bath of temperature T if the 
zero-point energy ^hv satisfies 



1 



-hv 



hv 



2 ew — 1 
This condition is equivalent to 

hv = TkT 



(14) 
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with r = In 3. For most of our considerations in the 
following, the value of the dimensionless constant T is 
irrelevant as the predictions are independent of it. 

Using eq. (fT5j) . the critical temperature T c in the 
Ginzburg-Landau model now corresponds to a critical 
frequency v c = TkT c /h where the gravitational activ- 
ity of photons ceases to exist. By putting eq. (TT5|) into 
eq. (|12p and (p~3|) we obtain our final result 



(16) 
(17) 
and 
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F{v) = -a -^ 1 - - 



valid for v < v c . For v > v c one has | V I / | 2 (^" 
F[v) = 0. 

Thus the number density |^| 2 of gravitationally active 
photons in the interval [v, v + dv] is nonzero for v < v c 
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only. In this way we obtain a finite dark energy density 
when integrating over all frequencies: 



where 



Pdark 



-hv\m 2 dv 
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(note the factor 1/2 as compared to previous work [IS . 
Il4j). The currently observed dark energy density in the 
universe of about 3.9 GeV/m 3 0, H || implies that the 
critical frequency v c is given by v c w 2.01 THz. 

Note that in our model virtual photons exist (in the 
usual quantum field theoretical sense) for both v < v c 
and v > v c , hence there is no change either to quantum 
electrodynamics (QED) nor to measurable QED effects 
such as the Casimir effect at high frequencies. The only 
thing that changes at v c is the gravitational behavior of 
virtual photons. This is a new physics effect at the in- 
terface between gravity and electromagnetism [ID, 
which solely describes the gravitational properties of vir- 
tual photons. 

We may calculate further interesting quantities for this 
electromagnetic dark energy model. The total number 
density N of gravitationally active photons is 



N = 



\^\ 2 {v)dv 
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Similarly, the total free energy density is given by 

f°° 327T 1 

Ftotai = J F{v)dv = ^^«o^. (23) 

Thus on average the free energy per gravitationally active 
photon is given by 



F 



total 
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This result is independent of v c and gives a simple phys- 
ical interpretation to the constant ao, which has the di- 
mension of an energy, just as for ordinary superconduc- 
tors. 

Our electromagnetic dark energy model depends on 
two a priori unknown parameters, the critical frequency 
hv c and the constant a . It shares many similarities 
with the Ginzburg-Landau theory of superconductors, 
formally replacing the number density of superconduct- 
ing electrons in the superconductor by the number den- 
sity of gravitationally active photons in the vacuum. It is 
instructive to see which values the constants ao and hv c 
take for typical superconductors in solid state physics. 

The Bardeen-Cooper-Schrieffer (BCS) theory yields 
the prediction [12j 
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Here /i denotes the Fermi energy of the material under 
consideration. For example, in copper /i = 7.0 eV, and 
the critical temperature of a YBCO (Yttrium-Barium- 
Copper Oxid) high-T c superconductor is around 90 K. 
This yields typical values of hv c ~ 8 • 10~ 3 eV and a ~ 
8 • 10-3. 

Remarkably, our dark energy model works well if the 
free parameters ao and hv c have the same order of magni- 
tude as in solid state physics. Many dark energy models 
suffer from the fact that one needs to input extremely 
fine-tuned or unnatural parameters. This is not the case 
for the Ginzburg-Landau-like model described here. Our 
model is based on analogies with superconductors, and 
in view of naturalness it would seem most plausible that 
the relevant dark energy parameters have a similar or- 
der of magnitude as in solid state physics. Moreover, the 
parameters of our model should be universal parameters 
related to electro- weak interactions, since we consider an 
electromagnetic model of dark energy. 

A possible choice is 
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where hv c ~ m v c 2 is proportional to a typical neutrino 
mass scale, and a e i « 1/137 is the fine structure con- 
stant. The motivation for (|27"|) is as follows. Since we are 
considering a model of dark energy based on electromag- 
netic vacuum energy, the relevant interaction strength 
should be the electric one described by a e i. Moreover, in 
solid state physics the critical temperature is essentially 
determined by the energy gap of the superconductor un- 
der consideration [TTI ] (i.e. the energy obtained when a 
Cooper pair forms out of two electrons). Something sim- 
ilar could be relevant for the vacuum. We could think 
that at low temperatures (frequencies) Cooper-pair like 
states can form in the vacuum. If this new physics has 
to do with neutrinos, one would expect that the rele- 
vant energy gap would be of the order of typical neutrino 
mass differences. Solar neutrino measurements provide 
evidence for a neutrino mass of about m v c 2 ~ 9 • 10 -3 eV 
[]~7L [l8j . assuming a mass hierarchy of neutrino flavors. 
This agrees with the energy scale hv c that we need here 
to reproduce the correct amount of dark energy density 
in the universe. 

Another constraint condition in our model of dark en- 
ergy is that the parameter a should not be too large. 
Otherwise, one would have in equilibrium a surplus of 
negative free energy density, which would counterbal- 
ance the positive dark energy density. We obtain from 
eq. ([22]), dill) and the ratio 



F 



total 



G4 



Pdark 105 

hence F tota/ | << p dark as required 



(28) 
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We now turn to possible measurable effects of our the- 
ory. The similarity with the Ginzburg-Landau theory 
of superconductivity, and in particular the fact that the 
potential parameters have the same order of magnitude, 
suggests the possibility that gravitationally active pho- 
tons could produce measurable effects in superconduct- 
ing devices via a possible synchronisation of the phases 
of the corresponding macroscopic wave functions. 

Denote the macroscopic wave function of gravitation- 
ally active photons by \&g (previously this was denoted 
as vp), and that of superconducting electrons (Cooper 
pairs) in a superconductor by \& s . So far we only dealt 
with absolute values of these wave functions, but we now 
introduce phases $g and & s by writing 



* s = |*»|e<*- (29) 
* G = |* G |e i<&G . (30) 

In superconductors one has \^ s \ 2 = \n S: where n s de- 
notes the number density of superconducting electrons. 
Similarly, in our model |\&g| 2 is proportional to the num- 
ber density of gravitationally active photons. Spatial gra- 
dients in the phase <I> S give rise to electric currents 
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where e is the electron charge and m the electron mass. 
Similarly, spatial gradients in the phase $g of gravita- 
tionally active photons would generate a current given 
by 

-* &fz tcfi 

JG = — |*g| 2 V* g = - — (*gV*g - 1-gWg)- (32) 
m Zm 

Whereas the strength of the electromagnetic current 
is proportional to the Bohr magneton /xg = -0-, the 
strength of the current given by eq. (|32|) is proportional 
to a kind of 'gravitational magneton' fia '■= ^ whose 
strength is a priori unknown. Presumably, /xg is very 
small so that this current is normally unobservable in 
the vacuum. 

In superconducting devices, however, the situation 
may be very different. Here both the phases <I> S and 
^G exist and the corresponding wave functions might 
interact. The strength of this interaction is a priori un- 
known since ^g represents new physics. If the interac- 
tion strength is sufficiently strong then in equilibrium the 
phases may synchronize: 



(33) 



This is plausible because of the similarity of the size of 
parameters of the corresponding Ginzburg-Landau po- 
tentials. If phase synchronization sets in, then fluctua- 
tions in <!>g would produce measurable stochastic electric 
currents of superconducting electrons given by 
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*s| 2 V$G- 



(34) 




v [THz] 



FIG. 1: Power spectrum S(y) := ^S(v) as given by eq. (|37[) 
at two different temperatures (solid and dashed lines) and 
comparison with the data of the experiment of Koch et al. 
[Tg| ] . Filled squares correspond to measurements at 4.2K, 
open squares to measurements at 1.6 K. 



However, these currents could only exist up to the criti- 
cal frequency v c . For v > v c one has $g — and hence 
V$g = 0. Noise currents that are produced by gravi- 
tationally active photons would thus cease to exist at a 
critical frequency v c given by about 2 THz. Generally, 
our Ginzburg-Landau model predicts that gravitationally 
active photons produce a quantum noise power spectrum 



S{v) = \hv ( 1 - V - 



(35) 



for v < v c and S(v) = for v > v c . 

In resistively shunted Josephson junctions, quantum 
noise power spectra induced by stochastically fluctuat- 
ing phases can be quite precisely measured [19] . For 
frequencies smaller than 0.5 THz the form of the power 
spectrum of current fluctuations has been experimentally 
confirmed [19j as 



a, x 4/1. hv 



(36) 



a direct consequence of the fluctuation dissipation the- 
orem [2(1 Hl[. Here R denotes the shunt resistor. The 
first term in eq. (|36|) is due to zero-point fluctuations 
and the second term is due to ordinary thermal noise. If 
there is full synchronization between the phases $ s and 
$G, then our Ginzburg-Landau model predicts a high- 
frequency modification of f|36[) given by 
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(37) 



This spectrum agrees with the spectrum (I36[) up to fre- 
quencies of about 1 THz but it then reaches a maximum 
at v max = fT x ^v c w 1.34 THz and approaches a cutoff 
at v c w 2.01 THz (see Fig. 1). New Josephson experi- 
ments are currently being carried out [22| that will reach 
the THz frequency range, thus being able to compare the 
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prediction (|3T[) with the experimental data for frequen- 
cies larger than 0.5 THz. 

To conclude, in this paper we have introduced a new 
model of dark energy where the dark energy of the uni- 
verse is identified as ordinary electromagnetic vacuum 
energy. The new physics of the model consists of a phase 
transition of virtual photons from a gravitationally active 
to a gravitationally inactive state. This phase transition 
is described by a Ginzburg-Landau-like model. The ad- 
vantage of our model is that it yields the correct amount 
of dark energy in the universe for natural choices of the 
parameters (quite similar to those used in the Ginzburg- 
Landau theory of superconductors), and that the predic- 
tions of our model can be tested by laboratory experi- 
ments. 

We end with a brief discussion of possible future de- 
velopments of the model studied here. First, it should 
be noted that our Ginzburg Landau approach can be 
applied to any type of boson and need not be confined 
simply to photons. Most generally we may assume that 
any particle in the standard model may exist in either 
a gravitationally active or inactive phase. To describe 
the phase transition behavior of a given boson, one may 
consider fundamental fermionic degrees of freedom that 
can condense into the boson being considered. A more 
advanced theory would construct the analogue of Cooper 
pairs in solid state physics, with a suitable weak attrac- 
tive force between the fundamental fermionic degrees of 
freedom. This would be equivalent to the development 
of a BCS type of theory that effectively reproduces the 
Ginzburg Landau model studied here. In 1 101 ] a simi- 



lar idea has been developed, with the weak attracting 
force leading to the condensate being gravity. In other 
models, for example for neutrino superfluidity, the weak 
mediating force between massive neutrinos that can lead 
to superfluid neutrino states is given by Higgs boson ex- 
change [23[ . Our model is phenomenological, just as in 
solid state physics the Ginzburg Landau theory arises as 
a phenomenological model out of the BCS theory. Our 
model does not rely on a particular microscopic model, 
but makes universal predictions for measurable spectra as 
shown in Fig. 1, assuming that the dark energy conden- 
sate interacts with ordinary Cooper pairs via a synchro- 
nization of phases. If other particle condensates also con- 
tribute to the measured spectra, then this would change 
the degrees of freedom entering into eq. ^ , thus leading 
to a slightly different critical frequency v c . Future exper- 
imental measurements of the critical frequency will thus 
be a very helpful tool to constrain the class of theoretical 
models considered. 
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